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We show that prominent features in voltage-current characteristics, recently measured in the mixed 
state of high-Tc superconductors and interpreted as evidence for an irreversibility line or a vortex- 
glass transition, may very well be explained with the simplest Kim- Anderson approach describing 
the vortex motion. In this case, the irreversibility line is not related to a transition in the system 
of vortices. Consulting numerous experimental reports on this subject we have not found a single 
example, which is in contradiction with this view. 

PACS numbers: 74.60.Ge, 74.60.Jg, 74.72.-h 



I. INTRODUCTION 

The commonly accepted picture of the magnetic flux 
motion in high temperature superconductors (HTSC) is 
extremely complex. Many different theoretical models 
have been developed in order to explain various aspects of 
the flux-creep process (see Ref. 1 for details). The com- 
plexity of the models that have been invented to capture 
the motion of vortices in HTSC is related to the apparent 
inability of the simple Kim- Anderson approach for ex- 
plaining the experimental results.Btl However, a serious 
reconsideration of available data indicates that this may 
not really be the case. Recent experiments have shown 
that the low-temperature flux-creep data may perfectly 
well be described using the Kim- Anderson approach, if 
a realistic .-profile of the pinning potential is taken into 
account J-!rla In this paper we use a similar approach to 
analyze the motion of vortices at temperatures close to 
the superconducting critical temperature T c . We show 
that the very specific features in voltage-current charac- 
teristics that manifest the flux-creep process and which 
are usually related to the irreversibility line and a vortex- 
glass transition, not only may be explained by employ- 
ing the Kim- Anderson approach, but are a direct conse- 
quence of this simple model. First, we briefly recapit- 
ulate the model and its consequences and subsequently 
address some issues concerning the irreversibility line and 
the vortex-glass transition. 
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FIG. 1. Schematic profiles of a pinning well for different 
values of the current density. The corresponding values of j/j c 
are indicated near the curves, (a) j/j c <C 1. Inset illustrates 
definitions of U\ and XJ?,. (b) 1 ~ j/jc <S 1. 



II. THE MODEL 



If the current density j in a sample is less than its 
critical value j c , all vortices are pinned and their motion 
is only due to a thermally activated hopping of the vortex 
lines over potential barriers or via quantum tunneling 
through the barriers. The probability of hopping in the 
case of thermal activation is 



w 



vq exp 



U 



(1) 



where v$ is the attempt frequency of the vortex line to 
surmount the potential barrier of height U and ks is the 



Boltzmann constant. 

As an example, a pinning-potential profile proportional 
to | sin x | is shown by the solid line in Fig. 1(a). At a given 
temperature T, the profile of the potential well for j = 
may be written as 



u(x) = U(T)f(x,T), 



(2) 



where U (T) is the energy fixing the pinning strength and 
the function f(x,T) defines the shape of the potential 
well, which may be temperature dependent as well. We 
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have chosen / such that \df/dx\ is maximum at x = 0, 
the site of the minimum of a single well. U(T) decreases 
with increasing temperature and it vanishes at the tem- 
perature at which the external magnetic field H is equal 
to the upper critical field H C 2- In the following, t/ max 
denotes the maximum value of u{x,j — 0) between two 
local minima as is illustrated in Fig. 1(a). Also U max is 
decreasing with increasing T but its temperature depen- 
dence may differ from that of U(T). 

An electric current does not change the interaction of 
a vortex with the pinning centers or other vortices, but 
it causes a Lorentz force Fl to act on the vortices. The 
Lorentz force tilts the potential profile, thus reducing the 
potential barriers in the direction of the vortex motion. 
In the presence of a current the potential profile may be 
written as 

u(x, j) = u(x, 0) — xFl (3) 

with Fl = j$o/c, j as the current density, <E>o the mag- 
netic flux quantum, and c the speed of light. The critical 
current density is reached if the potential barriers vanish. 
According to Eqs. (2) and (3), this results in 

«T) = ^2£, (4) 

where f' c is the maximum value of \df/dx\. As stated 
above, \df/dx\ reaches its maximum at x = 0, i.e., f c = 
\df /dx\ x=Q . Equation (4) is a formal definition of the 
critical current density, the only appropriate definition 
of j c that may be given in the mixed state of type-II 
superconductors. As will be shown in section IV (see Fig. 
6a), at temperatures close to the superconducting critical 
temperature T c , the sample resistance in the mixed state 
may be rather high even if j <C j c - This is the reason why 
different voltage criteria for the evaluation of the critical 
current density, although useful for practical purposes, 
are rather meaningless from the point of view of physics. 

The influence of current on the potential profile is il- 
lustrated in Figs. 1(a) and 1(b). It may be seen that 
at low currents (j <C j c ) the decrease of the activation 
energy with increasing current is entirely determined by 
the behavior of u(x) near its maximum [see Fig. 1(a)], 
while at currents close to j c only u(x) in the vicinity of 
x = is important [see Fig. 1(b)]. We denote the heights 
of potential barriers in the direction of the Lorentz force 
and opposite to it by U\ and U2, respectively [see inset 
to Fig. 1(a)]. For a non-zero current, U\ < TJi. The elec- 
trical field E in the sample is proportional to the average 
velocity of vortices and may be written as 

[ / Ui(T,j) \ ( U 2 (T,j) \] 

with 

S = -lv . (6) 



Here B is the magnetic induction in the sample and I is 
the vortex hopping distance. The second term in Eq. (5) 
describes the vortex hopping in the direction opposite 
to the Lorentz force. This term is, only important for 
3 /jc 1 and at high temperatures. 

If the function / in Eq. (2) is known, U% and U2 
may be deduced by using Eq. (3). In practice however, 
the situation is quite different. The main goal is rather 
to obtain information about u(x) from experimental re- 
sults. In most cases, the experimentally accessible infor- 
mation is limited to E(J) curves, measured at different 
temperatures. In order to obtain any useful information 
about u{x) from these experimental data, some a priori 
assumptions have to be made. As has been shown re- 
cently, the flux-creep rates in an epitaxial YF^CuaOy-z 
(YBCO) film may be very well described in a wide range 
of temperatures (0.02T C < T < 0.9T C ) with the assump- 
tion that not the shape, but only the amplitude of the 
u{x) function is temperature dependent In this case, 
the function / in Eq. (2) depends only on x and the tem- 
perature dependencies of t/ max and U (T) coincide. Based 
on this assumption, a scaling procedure .for the analysis 
of flux-creep data has been established.^ It is an impor- 
tant advantage of this approach that it allows for internal 
consistency checks and all assumptions that have been 
made may be verified retrospectively. It turns out that 
the scaling procedure breaks down at T > 80A" 0.9T c .eI 
This demonstrates that the assumption of a temperature 
independent / function is not valid close to the critical 
temperature and therefore should be abandoned in this 
temperature range. 

In the following we use Eq. (5) to analyze the voltage- 
current characteristics of the sample at temperatures 
close to T c . It will be shown that different consequences 
of the vortex motion, such as the vanishing of the per- 
sistent current at Ti rr < T c (irreversibility line) and 
the often observed sign change in the curvature of the 
log-E — log ,7 curves with increasing temperature, which 
is usually attributed to a vortex-glass transition, follow 
straightforwardly from Eq. (5) without any additional 
assumption. 

III. IRREVERSIBILITY LINE 

Among numerous unusual features of HTSC, the so- 
called irreversibility line (IRL) has attracted a lot of at- 
tention. This line in the H — T phase diagram separates 
two regions with distinctly different behavior Jalla Above 
the IRL the magnetization of the sample is reversible, 
which means that the sample cannot carry any persistent 
current. Below the IRL, irreversible magnetization is ob- 
served. Since the true superconducting state with zero 
dissipation is achieved only below the IRL, the position 
of this line in the H — T diagram is extremely impor- 
tant for applications. However, in spite of many years 
of intensive studies, the situation with regard to the IRL 
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is far from being clear. Not only the origin of the IRL 
is still under discussion, but there is also no consensus 
concerning valid procedures to establish the position of 
the IRL from the experimental data. 
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FIG. 2. Per- 

sistent current I p in a ring-shaped YBa2Cu307_ :E film as a 
function of temperature (description of the sample may be 
found in Ref. 4). The vertical arrow indicates the position of 
the irreversibility temperature as given by the present exper- 
iment. The sample was cooled in an external magnetic field 
H w 1000 Oe to T = 82 K. The current was subsequently 
induced by enhancing the magnetic field by 1 Oe. 



In order to make the physics more transparent, we con- 
sider here a ring-shaped sample. In this case the persis- 
tent current in the sample I p plays the role of an ir- 
reversible magnetic moment Mi rr . The current in the 
ring may easily be monitored by measuring the magnetic 
induction in the ring cavity using, for instance, a Hall 
probe placed in the center of the cavity. Typical exper- 
imental data for a ring-shaped sample obtained from an 
epitaxial YBCO film are presented in Fig. 2. This fig- 
ure shows a heating-cooling cycle of I p . Our data in Fig. 
2, as well as results of measurements of the irreversible 
magnetizationJjTJ reveal that above the irreversibility 
temperature, Ti rr , the persistent current is essentially 
zero. This is why Ti rr is usually considered as the temper- 
ature, at which the critical current density vanishes. It is 
commonly accepted that Ti rr corresponds to some tran- 
sition invoking the vortex system of the sample. Among 
a few considered possibilities for such a transition, the 
melting of a vortex-glass phase is the most .popular ex- 
planation for the irreversibility line in HTSCE3 However, 
as shown below, this kind of temperature dependence 
of the persistent current I p necessarily follows from the 
simplest Kim- Anderson approach for describing the flux- 
creep process. 

At temperatures well below the IRL the situation is 
quite clear. All vortices are pinned and their motion 
only occurs due to a thermally activated hopping of vor- 



tex lines over potential barriers. According to Eq. (1) the 
probability of hopping depends exponentially on V /ksT. 
Equation (1) implies that the resistance of a type-II su- 
perconductor in the mixed state formally never vanishes. 
However, at low temperatures and small currents, the ra- 
tio U /ksT is so large that the probability of hopping is 
negligible. In this case, an electric current induced in a 
ring-shaped sample may flow without a noticeable decay 
for years. Because the ratio lI/hsT decreases rapidly 
when approaching T c , the current-decay rate increases 
significantly with increasing temperature. 

The current decay rate dj/dt is proportional to the 
electrical field in the sample and therefore Eq. (5) may 
be used for evaluating this important quantity. In order 
to calculate the temperature and current dependencies of 
dj/dt, we have to assume some explicit expression for the 
profile of the potential well u(x). Unfortunately, practi- 
cally nothing is known about u{x) in this high tempera- 
ture range. In our previous study it was shown that the 
function / in Eq. (2) may be considered as temperature 
independent only for T < 0.9T c .eI For higher tempera- 
tures the temperature dependence of / has to be taken 
into account. In the following analysis we have chosen 
two rather different representations for f(x,T), i.e., 



f(x,T) = \x\-a(l-T/T c ) 



k x 2 



(7) 



and 

f{x,T) 

with 



x 



n I i \3/2 3/2 

(\x\ + x ) - V 



(1 - T/T c ) r 



U(T) = U {l-T/T c ) 3/2 



(8) 



(9) 



for both cases. This kind of U (T) follows from the sim- 
plest consideration of the vortex pinning at temperatures 
close to r c .E£EZl The potential profiles represented by Eqs. 
(7) and (8) are depicted in the insets to Fig. 3. We do 
not pretend that any of the / functions given by Eq. 
(7) or (8) exactly represent the real situation in HTSC. 
Our choice of / has been made in view of obtaining a 
reasonable profile shape and the possibility to perform 
analytical calculations. As will be shown below, the par- 
ticular choice of / does not influence the main qualitative 
features of the flux-creep process. 

We may now insert the potential profiles given by Eqs. 
(7) to (9) into Eq. (2) and calculate U x {T,j) and U 2 (T,j) 
with the help of Eq. (3). Substituting the resulting 
expressions for U± and TJi into Eq. (5), we obtain the 
corresponding temperature and current dependencies of 
the electrical field E in the sample. In order to com- 
pare the results of this calculation with the experimental 
data presented in Fig. 2, we need an evaluation of Eq 
and Uq entering Eqs. (5) and (9), respectively. For this 
purpose, we calibrated our calculated E{j) via a voltage- 
current characteristic measured for the same sample at 
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FIG. 3. Experimental E-j curve from Ref. 6 for T = 0.9T C 
and H — 0.9 kOe. The lines are approximations to this curve 
using Eqs. (7) and (8) for the profile of the potential well. 
The upper scale shows absolute values of thr current. The 
upper and lower insets show examples of the profiles of the 
potential well as given by Eqs. (7) and (8), respectively. 



J 



T = T w 0.9T c J2l Fig. 3 shows the experimental data 
together with the calculated E(j) curves using f(x,T) 
as given by Eqs. (7) and (8), respectively. In the first 
case, the constants Eq in and Uq were adjusted such as 
to approximate the low voltage part of the E(j) curve; 
the constant a = (1 — T /T c )~ k . In the second case, 
the parameters xq and b in Eq. (8) were used as fitting 
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FIG. 4. Calculated variation of dj/dt versus T/T c for 
f(x,T) given by Eq. (7) with k = 1.4. Calculations were 
made for 3 fixed current densities. The values of the renor- 
malized current densities j / j c (0.9T c ) are indicated near the 
curves. 
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FIG. 5. I p as a function of temperature. The solid lines are 
results of calculations for different representations of u(x). (a) 
On linear scales, the points represent the experimental data 
from Fig. 2. (b) logij, versus T/T c . 



parameters, as well. 

Using these adjusted expressions for E, we may now 
calculate different characteristics of the sample for the 
chosen f(x,T). Fig. 4 shows the temperature depen- 
dencies of dj/dt oc E for three values of the current den- 
sity. In this particular case, the calculation was made for 
f(x,T) given by Eq. (7) with k = 1.4. Note that we use 
a log-scale for the dj/dt-axis and that the total change 
in dj/dt is 100 orders of magnitude. This figure clearly 
demonstrates that dj/dt grows extremely fast with in- 
creasing temperature close to T c . This implies that a 
current which is practically constant in time at a certain 
temperature, will be decaying rather quickly at a slightly 
enhanced temperature. 

Experimentally the persistent current I p is usually de- 
termined as the current which is not decaying during the 
time of the experiment. This definition of the persis- 
tent current is not universal because the value of I p is 
obviously dependent on the experimental resolution and 
the time window of the particular experiment. In our 
approach we can calculate the temperature dependence 
of the persistent current by choosing dl/dt as set by the 
available resolution of the experiment. The following cal- 
culations were made for dl/dt = 3 • 10~ 4 A/s, which is 
the resolution for the experimental data presented in Fig. 
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2. In order to demonstrate how the I p (T) oc Mi rr (T) 
curve depends on the particular choice of the potential 
profile, the calculations were made for both approxima- 
tions of f(x, T) and for different values of the exponents 
k and to in Eqs. (7) and (8), respectively. The results of 
the calculations, together with experimental data from 
Fig. 2, are shown in Figs. 5(a) and 5(b). It may be 
seen that, depending on the degree of experimental reso- 
lution, an "irreversibility" temperature exists for all the 
chosen functions f(x,T). At the same time, both the 
shape of the I p (T) curves and the position of T irr are 
rather sensitive to the choice of /. This means that by 
a proper choice of U(T) and f{x,T) in Eq. (2), any 
experimentally observed temperature dependence of the 
irreversible magnetic moment Mi rr (T) or the persistent 
current I p {T) may be approximated to a high degree of 
accuracy and no specific transition in the vortex system 
is needed to explain the existence of the irreversibility 
line established in this way. 



IV. VORTEX-GLASS TRANSITION 

As has been mentioned above, the most popular in- 
terpretation of the irreversibility line in the H — T di- 
agram is ^fet-at is a manifestation of a vortex-glass 
transitionMSEJ Many experimental results seem to con- 
firm the concept of a voijtex-fflass melting at tempera- 
tures close to the IRL.II3E9ej In this section we show 
that the distinct variation of the shape of the logE— log j 
curves, which is usually interpreted as a manifestation of 
the vortex-glass melting, follows straightforwardly from 
our simple consideration outlined in the previous section. 

To make the point, we use f(x,T) as given by Eq. 
(8) with to = 2. The results of the calculations of E(j) 
are shown in Fig. 6(a) as E(j) curves at fixed temper- 
atures. Qualitatively this plot is indistinguishable from 
numerous experimental results (see Refs. 20-29) and it 
is clear that by corresponding adjustments of U(T) and 
/(#, T), any experimental E(J, T) curve may be approx- 
imated even quantitatively by this type of calculation. 
Fig. 6(b) shows the same E(j) curves as Fig. 6(a) at 
much lower voltages. It may be seen that the change of 
the curvature, which is usually attributed to the vortex- 
glass transition, is a universal feature of the E(J) curves 
at temperatures close to T c . With decreasing tempera- 
ture however, the sign change of the curvature is shifted 
to lower voltages where it is not accessible experimen- 
tally. 

We emphasize that the crossover to ohmic E(j) curves 
with decreasing current is independent of the particular 
choice of the u{x) function. As has been demonstrated 
in Section II, for j/j c <C 1, only a small part of the u[x) 
function near its maximum is important. In this region, 
u(x, 0) may be replaced by a Taylor expansion of the form 
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FIG. 6. (a) and (b) E—j curves calculated for f(x,T) given 
by Eq. (8) with m = 2 and parameters xq and b detirmined 
as is explained in the previous section. 



u{x) = U n 



2 V^o 



(10) 



Here xq is the point where u(x,0) has its maximum and 
A is the curvature of u at x = xq. Both A and xq may 
be temperature dependent. Using Eqs. (3) and (10), the 
activation energies are 



U U2 = U n 



it j 



. ^0*0 



f ( a:o$o 
2 \cAU ma , 



Ce^max 

Inserting (11) into Eq. (5), we get 

'{jxq^/cA) 2 



E = E exp 



exp 



2k B TU n 



sinh 



(11) 



\ ck B T 
(12) 



Because of the additional exponential term, the j depen-. 
dence of E cannot be reduced to a hyperbolic sinus .0 
Taking into account that Eq. (12) is only valid at the 
low-current limit, we may use the expansions of the ex- 
ponential and hyperbolic sinus functions and neglect all 
terms j n (n> 2). Eq. (12) is thus reduced to 



E 



j^-^E exp 
cksl 



Un 



kuT 



(13) 



and hence at any temperature, the sample resistivity p = 
E/j is independent of the current density for j <C Jo = 
cU ma , x /x $ . 
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A very similar explanation of the voltage-current char- 
acteristics near the " vortex- glassA transition has been 
suggested by Coppersmith at al.£3 In their short com- 
ment they considered a sinusoidal potential barrier. It 
was shown that even with this simple potential all the 
qualitative features of the experimental E(J) curves could 
be reproduced. The authors also pointed out that the 
insignificant quantitative disagreement with the experi- 
mental data is simply due to the arbitrary chosen sinu- 
soidal profile of the potential barriers. 



V. CONCLUSION 

As has been shown in Sections III and IV, the exper- 
imentally established vanishing of the persistent current 
at T = Ti rr < T C (H) does not necessarily mean that the 
critical current density vanishes at T > Ti rr and no sin- 
gular event, such as a vortex-glass transition, needs to be 
involved in order to explain the experimentally verifiable 
disappearance of the irreversible magnetization at the ir- 
reversibility line. This feature, as well as the change of 
curvature of the voltage-current characteristics with tem- 
perature, which is usually considered as a manifestation 
of the vortex-glass melting, automatically follow from the 
simplest consideration of thermally activated vortex mo- 
tion (see Figs. 5 and 6). 

In experimental papers dealing with vortex-glass melt- 
ing, there is usually no discussion of alternative possi- 
bilities to account for the experimental results. All the 
data are scaled according to the vortex-glass model and 
the corresponding correlation lengths are determined. A 
common statement in most of these papers is that the 
flux-creep model cannot provide a negatijae .curvature for 
the low temperature logE — log j curvesEj'CJ As may be 
seen in Figs. 6(a) and 6(b), a negative curvature may 
very well be accounted for by the flux creep model. In 
a treatment of thermally activated vortex hopping given 
in Section IV, this curvature depends on the particular 
pinning potential. 
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